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ABSTRACT 

We discuss the interpretation of distortions to stellar spectral lines with particular attention to line 
bisectors in the presence of an orbiting planetary companion. We present a simple model whereby light 
reflected by the companion can cause temporal variations to the observed line profiles. These distortions 
have a characteristic signature depending on the inclination angle of the system. For the known close-in 
extrasolar giant planets, the expected amplitude of the effect might not be far from current detection 
capabilities. This method could be used to directly detect the presence of the companion, yielding the 
orbital inclination and hence the planetary mass. Futhermore, a detection would measure a combination 
of the planetary radius and albedo, from which a minimum radius may be deduced. 

Subject headings: planetary systems — stars: individual (51 Peg, r Boo) — techniques: spectroscopic 
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1. INTRODUCTION 

The interpretation of low amplitude, periodic radial ve- 
locity variations in solar-type stars as being induced by 
orbiting low mass companions was recently put into ques- 
tion. Gray (1997) and Gray & Hatzes (1997) claimed to 
detect distortions to the line profile bisector (the locus of 
midpoints of a stellar absorption line from the core up to 
the continuum) in one spectral line (Fe i A6253) of the 
star 51 Pegasi at the now well known 4.23 day radial ve- 
locity period (Mayor & Queloz 1995; Marcy et al. 1997). 
They stated that the planet hypothesis could not account 
for variations in the line profiles, and proposed the alter- 
native explanation that 51 Peg was undergoing non-radial 
pulsations. Further observations and analyses (Gray 1998; 
Hatzes, Cochran & Bakker 1998a, b; Brown et al. 1998a, b) 
failed to confirm the claimed period in the line bisector, 
and the planet hypothesis has emerged as the most rea- 
sonable explanation of the radial velocity curve. 

Nevertheless, the question of interpretation of variations 
in the line bisector remains. Is it a given that the planet 
explanation is immediately excluded should intrinsic vari- 
ations in the line profiles of 51 Peg be detected at the 
claimed orbital period? Our purpose in this Letter is 
twofold; first, to provide a simple model whereby a close- 
in orbiting companion produces periodic distortions to the 
line profiles, and second, to explore the observational con- 
sequences of this effect, the amplitude of which might be 
close to current observational limits, and the specific signa- 
ture of which is set by the orbital inclination of the system 
and the radius and albedo of the companion. 

2. REFLECTED LIGHT 

The amount of light reflected by the planets of our solar 
system is a tiny fraction of the solar output. However, for 
the close-in extrasolar giant planets (CEGPs) discovered 
to date, this fraction increases substantially due to the 
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proximity of a large planet to the star. The semi-major 
axes of these systems, r Boo, 51 Peg, v And and p 1 Cnc, 
are 0.046, 0.051, 0.057 and 0.11 AU respectively (Mayor 
& Queloz 1995; Butler et al. 1997). An estimate of the 
flux ratio can be achieved as follows. Let D denote the 
distance from the star to the planet, R p the radius of the 
planet, R s the radius of the star, and a the angle between 
the star and the Earth as seen from the planet. We assume 
that R p <C R s <C D and that each surface element on the 
planet reflects a fraction rj of the incident flux uniformly 
back into the local half-sky (Lambert's Law). If r\ does 
not vary with position, then it is simply the Bond albedo. 
Integration over the sphere yields the ratio of the observed 
flux of the planet at a = to that of the star, 



(1) 



e= iU 

For either 51 Peg b (D = 0.051 AU and assuming R p ~ 
1.3i?j up (Guillot et al. 1996)) or r Boo b (D = 0.046 AU 
and assuming R p ~ 1.2i?j up ), the result is e ~ r\ x 10~ 4 . 

The observed flux ratio will vary with the angle a, which 
is a combination of the orbital inclination i and the orbital 
phase cf>. Here we define the orbital phase to be the one 
conventionally given by radial velocity measurements; for 
a circular orbit (the only case we consider, as tidal effects 
circularize the orbits of these extremely close companions), 
<p = is taken to be the time of maximum recessional ve- 
locity of the star. The angle a € [0, ir] is then given by 
cos a = — sini sin27r0. Integration of the intensity over 
the surface of the sphere viewed at an angle a yields the 
phase dependent observed flux ratio, 



sin a + (it — a) cos a 



(2) 



For a discussion of the albedos and phase functions of so- 
lar system objects, see Harris (1961). Note that in the 
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expression for /(</>, i), we have ignored the possibility of 
occupation; eclipses would have a much more pronounced 
effect and are ruled out for r Boo, 51 Peg and p 1 Cnc 
(Henry et al. 1997; Baliunas et al. 1997). 

Thus CEGPs at moderate to high inclination and of high 
reflectivity will produce a photometric modulation of the 
system which will be accessible to upcoming satellite mis- 
sions, as discussed by Charbonneau (1998). However, pho- 
tometry does not make use of the spectroscopic separation 
of the planet and star. The key to the method we propose 
here is in analyzing the spectrum of the system. Assum- 
ing the planet reflects an essentially undistorted spectrum 
of the star, a copy exists in the observed stellar spectrum 
with relative amplitude /(</>, i) but at a location which is 
Doppler shifted due to the motion of the planet. 

The ratio of the amplitude of the radial velocity varia- 
tions of the planet to that of the star, K p /K s , is simply 
the ratio of the mass of the star to the mass of the planet, 
M s /M p . From the single line spectroscopic orbit, only 
the minimum mass of the planet, M p sini, is determined 
(assuming the stellar mass is estimated via other means). 
Thus, only the combination K p / smi is known. For 51 
Peg, K s ~ 56 ms" 1 and M p smi ~ 4.4 x 10" 4 M S , so 
K p J sinz ~ 130 kms -1 . Thus there can be a huge differ- 
ence in velocity amplitude between the stellar and reflected 
light signal. 

For a given stellar spectral line, the reflected line will 
move back and forth as the phase changes. At <j> — 0.25 
and (f) — 0-75 the radial velocities of the star and planet 
will be the same and the reflected line will lie on top of the 
stellar line. At cj) = and <f> = 0.5, the lines will be maxi- 
mally separated. A single spectrum with sufficient signal- 
to-noise would in principle be able to show both. We have 
obtained many high signal-to-noise, high resolution spec- 
tra of t Boo, and are attempting to isolate the Doppler 
shifted reflected spectrum from the orders of magnitude 
brighter stellar spectrum by making use of the known or- 
bital parameters supplied by radial velocity measurements 
(Butler et al. 1997). The results of this complementary in- 
vestigation will be presented in a separate paper. Here we 
present an alternate method of detection which capitalizes 
on the high resolution observations and analysis methods 
that have already been developed for studying line bisec- 
tors. The motion of the reflected line through the stellar 
line will cause a small distortion which might be detectable 
by careful measurements of the line bisectors or through 
other techniques. 

3. MODEL AND RESULTS 

We assume that the planet reflects an undistorted spec- 
trum of the primary, leading to a composite spectrum con- 
sisting of a reflected component greatly reduced in am- 
plitude and Doppler shifted by the orbital motion of the 
companion. We model this process for an individual iso- 
lated spectral line. Typical line bisectors for solar type 
stars show a distinctive C shape caused by granulation 
in the photosphere (see e.g., Gray 1992). Here, we have 
described the lines as simple Gaussian profiles of appro- 
priate width and depth for the star in question. While we 
have calculated more detailed models (multiple Gaussian 
components to the line) so as to reproduce the constant 
bisector C shape, the effect of this is simply to add a time 
invariant term to the bisector shape derived from the sym- 



metric Gaussian model. 

The velocity of the planet relative to the star at phase 
<j) is simply 



v p (<j),i) = -K, 



M s + M p 
cos 2wo 



(3) 



We have modelled the stellar and planetary line profiles as 
a function of velocity by 

s(v) = l-a s exp(-v 2 /Av s 2 ) (4) 
p{v,(j>,i) = 1 - a p exp(-(w - v p {(j), i)) 2 / Av p 2 ) (5) 

where we have explicitly allowed for a different planetary 
line width {Av p ) and depth (a p ) due to the rotation rates 
of the planet and star (see below) . 

The flux ratio is given by Equation ^, thus the observed 
spectrum of the system is 



I(v,ct>,i) 



s{v) + f{<j>,i)p{v,<l>,i) 
1 + 



(0) 



where we have renormalized to the continuum flux level. 
The effect of the reflected light component on the observed 
bisector is shown schematically in Figure [l]. 

For each point in the phase, we deduce the bisector by 
evaluating the midpoint between the two halves of the ab- 
sorption feature at a number of flux levels, interpolating 
between sampled points to the desired flux. We have also 
characterized these distortions in terms of two parameters 
conventionally used in bisector analyses, the velocity span 
and the bisector curvature. The velocity span measures 
the velocity difference in the bisector at two flux levels; 
we have chosen these levels at two percent below the con- 
tinuum (/ = 0.98) and two percent above the line depth 
(/ = 1.02 — a s ). The bisector curvature is defined to be 
the velocity span of the top portion (between / = 0.98 
and 1=1 — (a s /2)) minus the velocity span of the bottom 
portion (between 1=1 — {a s /2) and / = 1.02 — a s ) of the 
bisector. 

For convenience, we assume r\ = 1 in our model calcu- 
lations. The amplitude of the distortion to the bisector 
(and hence to the velocity span and curvature) is linear in 
r\. The results we present below should therefore be scaled 
for a given choice of the reflectivity. 

We have computed simulations for 51 Peg with e = 10~ 4 
and a stellar line width of Av s = 2.2 kms -1 . We as- 
sume the reflected light is undistorted, so that Av p = Av s 
and a p = a s . We sample the data simulating a spectral 
resolution A/AA = 200,000 comparable to current high- 
resolution observations. The predicted distortions to the 
bisector for i = 80° are shown in the left panel of Figure 
|2| for a selection of radial velocity phases. The deduced 
velocity span and bisector curvature are shown in the top 
panels of Figures 3 and 4 respectively. The bisector is 
distorted for the fraction of the phase where the stellar 
and planetary lines overlap and the planet is sufficiently 
illuminated as viewed by the observer. For a near edge-on 
geometry in the case of 51 Peg, this corresponds to roughly 
3% of its 4.23 day period, or about 3 hours, centered on 
4> = 0.75. Although the stellar and planetary lines also 
overlap at <fi = 0.25, there is less distortion as the planet 
is less illuminated as viewed by the observer. For lower 
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inclinations, the fraction of the period in which there is a 
significant distortion is greater, but the amplitude of the 
distortion is less. 

We have also modelled this effect for another CEGP sys- 
tem, r Boo (K s = 469 ms -1 and M p sini/M s = 3.1 x 1(T 3 , 
(Butler et al. 1997)). Observations of Ca n H + K lines 
suggest that the rotational period of the star is the same 
as the observed orbital period of the companion (Baliu- 
nas et al. 1997), implying that the star and planet are 
tidally locked. The synchronization timescales for r Boo 
are discussed in Marcy et al. (1997). If indeed the system 
is tidally locked, the planet does not see the rotational 
broadening of the spectral lines of the primary, as there 
is no relative motion between any point on the surface 
of the planet and any point on the surface of the star. 
Thus the planet reflects a rotationally unbroadened spec- 
trum of the primary, with an intrinsic width due to the 
stellar photosphcric convective motions. An observer then 
sees this line broadened only by the planetary rotation, 
which in this tidally locked scenario is small. Thus in the 
case of r Boo, we predict relatively narrow planetary lines 
(we adopt Av p ~ 4 kms -1 ) superimposed on the much 
broader stellar lines (Ad s ~ 15 kins" 1 ). The results of 
our simulated bisector distortions for i = 80° are shown 
in the right panel of Figure || for a selection of radial ve- 
locity phases. The deduced velocity span and curvature 
with phase are shown in the bottom panels of Figures |^ 
and [| respectively. The line is significantly distorted for 
the fraction of the phase where the planetary and stellar 
lines overlap, which now is substantially increased due to 
the high degree of rotational broadening experienced by 
the lines of the star. For an edge-on configuration, this 
amounts to roughly 10%, or about 8 hours of the 3.3 day 
period, centered on cj) = 0.75. Due to the different line 
widths, the deduced bisector varies with a much greater 
amplitude and in a more complex manner than in the case 
of 51 Peg. The velocity span and curvature do not well 
characterize the behavior of the bisector. We advocate 
that the full bisectors be examined in a search for this 
effect. 

4. DISCUSSION 

We have shown that the presence of a close-in extrasolar 
giant planet can distort the observed spectral line profiles 
of the system through the contribution of reflected light 
from the companion. For 51 Peg, assuming i — 80° and 
e = r\ x 10~ 4 , the amplitude of the distortion is roughly 
Arj ms -1 and 3r/ ms" 1 (peak-to-peak) to the velocity span 
and curvature, respectively. If rj is high, the effect is com- 
mensurate with current observational limits: Hatzes et al. 
(1998a) combined measurements of five spectral lines of 
51 Peg and achieved a scatter of 7 ms -1 and 12 ms -1 
for the velocity span and curvature. More recent results 
(Hatzes et al. 1998b) achieve a precision of 1 ms -1 for the 
velocity span and 4 ms -1 for the curvature. In the case 
of t Boo, assuming that the system is tidally locked and 
that i = 80° and e — rj x 10 -4 , the velocity span and cur- 
vature distortions are both about 42ry ms -1 peak-to-peak. 
Again, for a high value of rj, these are close to the current 
upper limits as measured by Hatzes & Cochran (1998). 

The amplitude of the velocity span and curvature vari- 
ations we have stated here rely upon measurement of the 
line bisector to within 2% of the continuum flux level. Bi- 



sector measurements tend to avoid the extreme wings of 
the line profile as dl/dX is small and hence the deduced 
bisector is more easily corrupted by noise. For 51 Peg, 
the equal lincwidths of the stellar and reflected spectrum 
mean that much of the signal in the velocity span and 
curvature is a result of the distortion of the bisector close 
to the continuum, as shown in Figure 2. For r Boo, the 
different line widths in the tidally locked scenario produce 
a second trough and peak in the velocity span and cur- 
vature which are insensitive to the height to which the 
bisector is measured, because the distortion occurs away 
from the continuum. Note that confirmation of the pre- 
dicted distortions does not require accurate measurement 
of the intrinsic line bisector, only precise measurement of 
the temporal variations to it, and thus it is not threatened 
by contamination from line blends. Also, isolated lines 
are not required as distortions induced by the transit of 
nearby lines can easily be modelled under the assumption 
that the planet reflects a Doppler shifted copy of the stellar 
spectrum. 

The estimate of the reflected light amplitude presented 
in Equation 1 does not consider the underlying source of 
the albedo, nor any angle dependent scattering effects as- 
sociated with it. Estimates of the albedos of these CEGPs 
are uncertain as the chemistry of these atmospheres is un- 
known. The gas giants of the solar system have Bond 
albedos ranging from 0.7 to 0.9. The large reflectivity is 
primarily due the presence of ices and clouds high in the 
atmosphere, and these could not exist at the much greater 
effective temperatures of CEGPs. Dominant sources of re- 
flectivity for CEGPs are scattering from dust and Rayleigh 
scattering. Detailed modelling of the atmospheric chem- 
istry of these planets indicates that condensates such as 
MgSiOa may indeed form at the required height (Burrows 
et al. 1997). Recent models by Seager & Sasselov (1998), 
which explicitly solve the radiative transfer through a 
model atmosphere for t Boo b, predict an albedo of less 
than 1% at 5000 A. Although this model is incomplete in 
that it considers only one dust species, of interstellar size 
distribution and solar adundance, it seems unlikely that 
albedos as high as those of the solar system gas giants 
could be reached by a dusty CEGP; further investigation 
is called for. An upper limit on the bisector distortions 
which is sufficiently strict to disallow a large value of rj 
would provide support for such models. 

5. CONCLUSION 

Via reflected light, a close-in extrasolar giant planet 
could produce a distortion to the stellar spectral line bi- 
sectors at a level which is approaching the current obser- 
vational upper limits, given a sufficiently high albedo. The 
variation of the bisector distortions with phase has a very 
distinctive signature which depends on the inclination of 
the orbit of the planet to the line of sight. Measurement of 
such periodic bisector variations would constitute a direct 
detection of the companion, a major advance in the field 
of extrasolar planet research. This would yield the incli- 
nation, and hence the mass of the companion. It would 
also measure the product of the albedo and the square of 
the radius of the planet, from which a minimum planetary 
radius may be deduced. An estimated albedo would yield 
a value for the planetary radius, which can be combined 
with the measured mass to calculate such critically impor- 
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tant quantities as the surface gravity and average density. 
Conversely, a firm upper limit on such distortions would 
place a constraint on the albedos of these planets. Either 
by detection or by an upper limit, an accurate measure- 
ment of the line profile bisector would place an important 
observational constraint on models of these close-in extra- 
solar giant planets. 
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Fig. 1. — The observed line profile is the sum of the stellar absorption feature and a time varying, Doppler shifted reflected light component. 
The amplitude of the reflected light has been greatly exaggerated in this plot so as to make the resulting distortion to the bisector visible. 
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Fig. 2. — Predicted distortions to the bisector of a typical spectral line of 51 Peg and r Boo for i = 80° and e = 10 4 . The bisectors are 
labelled by radial velocity phase. 
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Fig. 3. — The velocity span variations with radial velocity phase for 51 Peg and t Boo, both with e = 10 4 . The inclinations are 80° 
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Fig. 4. — The curvature variations with radial velocity phase for 51 Peg and t Boo, both with e = 10 4 . The inclinations are 80° (solid), 
25° (dashed) and 10° (dotted). 



